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ABSTRACT 

The excitation of solar and solar-like g modes in non-relativistic stars by arbitrary 
external gravitational wave fields is studied starting from the full field equations of 
general relativity. We develop a formalism that yields the mean-square amplitudes and 
surface velocities of global normal modes excited in such a way. The isotropic elastic 
sphere model of a star is adopted to demonstrate this formalism and for calculative 
simplicity. It is shown that gravitational waves solely couple to quadrupolar spheroidal 
eigenmodes and that normal modes are only sensitive to the spherical component of 
the gravitational waves having the same azimuthal order. The mean-square amplitudes 
in case of stationary external gravitational waves are given by a simple expression, a 
product of a factor depending on the resonant properties of the star and the power 
spectral density of the gravitational waves' spherical accelerations. Both mean-square 
amplitudes and surface velocities show a characteristic i?*-dependence (effective R^- 
dependence) on the radius of the star. This finding increases the relevance of this 
excitation mechanism in case of stars larger than the Sun. 
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1 INTRODUCTION 



Solar-like oscillations in non-relativistic stars are thought 
to be stochastically excited by turbulent convection 
in the o uter layers of the star. Theoretical inves- 
tisations (iGoldreich fc KeelevI ITqTtI: iBalmforthl 1 19921 : 



' iGoldrcic h. Murray fc Kumailll994l : ISamadi fc Goupilll200ll ') 
identified two major mechanisms by which oscillations are 
driven: the first source term is due to the turbulent Reynolds 
stress and is thus related to the Reynolds stress tensor; the 
second mechanism stems from the advection of Eulerian en- 
tropy fiuctuations by turbulent m otions (the en tropy source 
term). A generahzed formalism bv lSamadi fc Go upil (2001) 
for radial modes, which attempts at combining both mech- 
anisms in a consistent force d wave equation, was further ex- 
tended by iBelkacem et al] (2008) to the case of non-radial 
modes. Theoretically computed excitation rates for radial 
solar p modes are in good agreement with observational 
data f rom the GOLF instrument on board the SOHO space- 
craft (|Samadi et al.ll2003al lbl: [Belkacem et al.ll2006al lbl). Fur- 
thermore, radial solar p- mode peak heights were s hown to 
be re asonably well modelled jChaplin et al.l l2005l : iHoudelj 
I2OO6I ). However, in case of solar g modes the entropy source 
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term is negligible (|Belkacem et al.l |2009| ) and quantitative 
estimates of mode amplitudes differ from each other by 
orders of magnitude, depending mainly on the assumed 
eddy-time correlation functio n, that is on how the turbu- 
lent eddies arc time-correlated (|Kumar. Quataert fc Bahcalll 
[l996: Belkaccm et al. 20091: see in particular the discussion 
in I Appourchaux et al.ll2O10l ') . Since the detection of solar g 
modes still remains controversial, theoretical predictions of 
their amplitudes are of utmost importance. 

In the present paper, we consider an independent driv- 
ing mechanism which we expect to be relevant for solar and 
solar-like g modes, especially in case of stars larger than the 
Sun that are nearby strong gravitational wave sources: ex- 
citation by external gravitational waves. Due to damping 
rates, I = 1 and I — 2 so lar g modes are the m ost probable 
candidates for detection (|Belkacem et al.ll2009l ). As is shown 
below, this mechanism is able to excite quadrupolar eigen- 
modes {I = 2) and thus it may infiuence the future detection 
of solar and solar-like g modes. 

The excitation of normal modes by gravitational 
waves has already been studied for relativistic stars. 
These stars were shown to have a family of normal 
modes that are directly associated with the curvature of 
space-time, called gravitational- wave modes (w modes), 
and th at have no Newton i an an al ogue in non-relativistic 
stars l|Kokkotas fc Schut j 1 19921 ). lAndersson fc Kokkotaj 
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l|l996l ) first showed that these normal modes can be 
excited by gravitational waves impinging on the star. 
This initial work was extended and refined by many 



other studies (see e.g. the review b y Kokkotas fc Schmidtl 
1 19991 . iBernuzzi. Nagar fc De Pietril |2008| and references 
cited therein). Moreover, since the mid-1990s astero- 
seismology and gravitational wave physics were com- 
bined in order to gain inf orm ation a bout the interior 
of relativistic stars (see e.g . lAndersson fc Kokkotas! 1 19981 : 
iBenhar. Ferrari fc Gualtieril |2004| ). Even in case of non- 
relativistic stars such as the Sun, gravitational wave signals 
generated by g modes have been studied hoping that a future 
gravitational wave detector such as the Laser Interferometer 
Space Antenn a (LISA) will be able to unequivocally detect 
solar g modes (|Polnarev. Roxburgh fc Baskaranll2009l ). 

However, a thorough theoretical investigation of an ex- 
citation mechanism for normal modes of non-relativistic 
stars due to gravitational waves has not been presented so 
far and is th erefore worth studying even for its own theoret- 
ical interest. iKhosroshahi fc Soboutil (119971 ) already consid- 
ered this excitation mechanism, but concentrated on calcu- 
lating energy absorption cross-sections. Other authors stud- 
ied such a driving mechanism focusing on the excitation of 
elastic waves in solid planets or on the excitation of nor- 
mal modes of resonant spheres used in gravitational wave 
detectors. In the present paper, however, we take a differ- 
ent perspective: motivated by the aim of calculating surface 
velocities for stars, which are interesting quantities for aster- 
oseismology, we present a new way of solving this excitation 
problem. 

The general idea of the present paper is to combine the 
physics of sp herical resonant mass gravitational wave detec- 
tors (see e.g. lAshbv fc Dreitlein|[l975l : lLobdll995l : iMaggiord 
|2008| . and references cited therein) with asteroseismology, 
i.e. we consider a star as a giant gravitational wave detec- 
tor. The equation of motion for the internal velocity field of 
a star under the action of arbitrary external gravitational 
waves is derived and solved analytically. This analytical so- 
lution to the equation of motion enables us then to derive 
analytical expressions for the surface velocities, which in 
turn are crucial for assessing the astrophysical and aster- 
oseismic relevance of the proposed excitation mechanism. 
However, we do not conduct a hydrodynamical calculation, 
but instead we take an isotropic elastic sphere as a very 
simple and crude model of a star for calculative simplicity. 
This resonant sphere model is ideally suited to demonstrate 
the formalism of the proposed excitation mechanism and it 
eventually yields the same analytical results as in a hydro- 
dynamical computation. This is due to the fact that there is 
a similar equation of motion in both cases and it is also due 
to the spherical symmetry of the problem, since, as it turns 
out, the separation of the star's spheroidal normal modes 
into a horizontal and a radial component together with the 
orthogonality relation of the normal modes already consti- 
tutes the relevant information about the normal modes that 
is needed to analytically solve the equation of motion for the 
internal velocity field. Employing the elastic sphere model 
also makes this formalism applicable to the excitation of 
elastic waves in solid planets, which, however, is not the fo- 
cus of the present paper. 

The paper is organized as follows. In Section 2 we 
briefiy discuss the adopted approximation of linearized grav- 



ity and some properties of gravitational waves in transverse- 
traceless gau ge that are needed for th e subsequent analy- 
sis. Following^ Ashbv fc DreitleinI (| 19751 ). we start from Ein- 
stein's equations of general relativity and the Bianchi identi- 
ties in order to derive the inhomogeneous equation of motion 
for the internal velocity field of the star under the action of 
arbitrary external gravitational waves in Section 3. In this 
section we also discuss in detail the various assumptions and 
approximations that enter the computation. Section 4 is de- 
voted to a very brief discussion of the normal modes, which 
are the solutions to the homogeneous equation of motion 
without damping and forcing terms. In Section 5 we analyt- 
ically solve the inhomogeneous equation of motion and cal- 
culate the square amplitudes of the excited normal modes. 
These square amplitudes finally allow for an analytical com- 
putation of surface velocities in Section 6. Conclusions and 
prospects are presented in Section 7; an appendix is added 
to shorten the calculation in Section 5. 



2 GRAVITATIONAL WAVES 

For the present analysis we assume that the linearized theory 
of general relativity is applicable, so that the metric takes 
the form 

where \h^,^\ -C 1 and 77^,^ = diag(— 1, -|-1, -|-1, -1-1). As is 
shown in Section 3, we are free to choose the transverse- 
traceless (TT) gauge, defined bjQ 

h""" = 0, h\ = 0, ^'h^J = 0. (1) 

The gravitational waves incident on the star under consid- 
eration satisfy the linearized field equations in vacuo, which 
in TT gauge read 

Uhij = 0. 

From equation ([T}, it is obvious that hij is a spin- two 
tensor field and thus in an expansion in spherical harmonics, 
there are solely contributions from I = 2: 



2 

E 

m=-2 



(2) 



where hm are the five independent spherical components 
of hij, and X^/* Bxe components of the following matrices, 
which form a basis for the five-dimensional space of traceless 
symmetric tensors: 
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V^" = W — — 



^ Throughout this paper, Greek indices take space-time values 
0,1,2,3, whereas Latin indices take spatial values 1,2,3 only. Re- 



peated indices are summed over. 
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These matrices are linked to spherical harmonics via 

Y2m{0,(j}) = yfj^UtUj, 



(3) 



where Ui are the components of the radial unit vector n = 
(sin 6* cos sin ^ sin 0, cos and the orthogonality relation 
is given by 



Stt 



(4) 



(see, e.g., lMaggiord[2003 ). 

For later use we mention two additional important re- 
lations. Multiplying equation by niUj, summing over i 
and j and inserting equation ((3|, we obtain 



m=-2 



hmY2m{6, 4>), 



(5) 



and inverting equation ((3]) with the help of equation Q 
yields 



"i'T-i ~ ^Sij = ^ C™ Y2m.{ 
m — — 2 



(6) 



where the factor of 1/3 is fixed by the requirement that the 
left hand side be traceless, since the coefficients c^j are given 



3 EQUATION OF MOTION 

In order to derive the equation of motion for the internal ve- 
locity field of a star under the influence of external gravita- 
tional waves, we start with the full field equations of general 
relativity and the Bianchi identities, which together imply 
the conservation equations of energy and momentum. 



0, 



(7) 



where T''" are the components of the stress-energy tensor 
of the star under consideration. The equation of motion is 
obtained from the spatial components of equation 0: 



+ 



J- MP-' 



*- MP-' 



Since the centre of mass of the star will move on a 
geodesic in space-time, it proves useful to work in Fermi 
normal coordinates with the origin at the centre of mass at 
all times. In this reference frame 

rf,,[p(t)] = 0, 9or^.[PWl = o, 

where Pit) denotes the centre of mass at time t. Using these 
expressions one finds for the Riemann tensor 

Alternatively, within the framework of the linearized 
theory - which we assume is valid here - the Riemann ten- 
sor is invariant, rather than just covariant, and thus can be 
evaluated in any preferred frame. Consequently, choosing 
the TT frame for convenience one has from the linearized 
theory the following expression for the components i?'ojo of 
the Riemann tensor in terms of the metric: 



where hij are the gravitational wave components of the met- 
ric in TT gauge. 

It is now assumed that the diameter d of the star is much 
smaller than the typical length scale A/27r over which the in- 
cident gravitational radiation changes substantially. Clearly, 
this is an assumption that in practice has to be checked case 
by case. Due to the fact that under this assumption the com- 
ponents hij have essentially no spatial dependence over the 
volume of the star, we obtain the following relation 

which can be integrated to give 



r nn — 



2c2 



Moreover, we assume that the internal motions of the 
star are non-relativistic. In this Newtonian approximation, 
only r"" terms need be retained on the right-hand side of 
equation © and T''" is given by T^" = pc^, where p is 
the equilibrium proper mass density of the star. Identifying 
T^^ /c as the non-relativistic momentum density given by 
pv, and T*-* as the negative non-relativistic stress tensor. 



we arrive at 



dvi do 



dt 



(9) 



where v with components Vi denotes the internal velocity 
field of the star. 

For the reasons stated above in Section 1, the star is 
treated as an isotropic elastic sphere. This requires p{x) — 
p{r). Due to the external force exerted by the gravitational 
waves, an infinitesimal volume element of the elastic sphere 
centred at position x will be displaced according to cc -|- 
u{x, t), where we assume the displacements to be sufficiently 
small such that the linear theory of elasticity is appropriate. 
To this approximation and neglecting self-stresses caused 
by the intrinsic gravitational field (see also the discussion 
in Section 4), the elastic stress tensor for isotropic media is 
given by 



(10) 



where uim = il/2)( diUm + dmUi) and A and p are the usual 
Lame coefficients jLandau fc LifshTt3 Il963l ). The positive 
constants A' and p' parametrize the viscous properties of 
the elastic sphere. Note that equation (jlOp only holds in 
Fermi normal coordinates. 

Since the velocity field of a star is the easier measurable 
quantity than the displacements themselves, e.g. by Doppler 
techniques, we differentiate equation ® with respect to time 
in order to obtain the equation of motion in terms of the 
velocity field: 



C v + V{v) 



df 

dt ' 



where f{x,t), with components given by 
1 

2' 



-phijx\ 



(11) 



(12) 



denotes the external driving force per unit volume exerted 
by the gravitational waves. The linear differential operators 
£ and D act on arbitrary differentiable vector fields X ac- 
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cording to 

C[X) = i f(A + ^)V(V-X) + MV^A:l 

p '- 

V{X) = -(A' + ^')V(V-X)-^t'V^X. 

The model of viscosity adopted here corresponds to 
the Kelvin- VoiRt mod el considered by ^Ortega (1992.) and 
iLobo &: Ortegal l|l998l ). For the present paper, however, it 
suffices to consider the damping operator T) as some general 
abstract operator that is linear in its argument and does not 
involve any time derivatives. These are the only properties 
that enter and affect the calculation conducted in Section 
5. Hence, in principle, any other (phenomenological) model 
may be employed, provided that the resulting damping op- 
erator has the aforementioned properties. 



4 NORMAL MODES OF AN ELASTIC 
SPHERE 

The normal modes of an elastic sphere form a basis of so- 
lutions to the homogeneous equation of motion without a 
damping term, 



£1 

dt2 



0. 



(13) 



This equation of motion has to be supplemented with the 
appropriate boundary condition 







(14) 



at r — R, where A' and n' are set to zero, = Xi/r is 
the unit normal in Cartesian coordinates and R denotes the 
radius of the sphere. In order to find the normal modes, 
one can assume the displacement to have a harmonical time 
dependence, 

it(a;,t) = *(a;)e-'"*. (15) 

Inserting equation psp into equation p3p yields the eigen- 
value problem 



(16) 



Solving for the eigenfunctions and eigenfrequencies of £ sub- 
ject to the boundary conditions given by equation (|14[) is a 
standard problem in elasticity theory und thus shall not be 
reconsidered here. We solely quote the final result below for 
later use in Section 5 (cf. e.g. lAshbv fc Dreitleinlll975l : iLobd 
ll995l : lMaggi"or3l2008 h. 

There are two families of solutions to equation (|16|l . The 
spheroidal modes can be written as 



* , (r 9 



■ a„i(r)Ytmer - bni{r)ier xL Yi„ 



(17) 



where L = —ix X V is the angular momentum operator. 
The normalization is fixed by requiring 



d'a;p(*f,„) 



(18) 



where I is normally set equal to the mass M of the sphere 
and V denotes its volume. For I — 0, the spheroidal modes 
are purely radial, that is, b„oir) = 0. The toroidal modes 
are given by the expression 



where I > 1, and are normalized in the same way as 
the spheroidal modes. The orthogonality and normalization 
properties of the eigenfunctions of £ can be summarized as 



IS. 



(20) 



where = {nlm; SoiT} is an abridged index. It is impor- 
tant to stress that in the present paper we are not interested 
in detailed analytical expressions for the coefficient functions 
a,ni{r), bni{r) and c„i(r), we rather treat them as variables, 
i.e. we express the final results as functions of these coeffi- 
cient functions. 

In deriving equation (|10[) . effects of self-stress caused 
by the intrinsic gravitational field of the sphere were ne- 
glected. However, these effects are expected to be non- 
negligible for spher es as massive as astronomical objects 
(.Ashbv fc DreitleinilTgTSi '). Including these effects into the 
present formalism yields a different differential operator jC 
in equation (|11|) . 

£(X) = — {{X + ■ X) + fiV^X 

po ^ 

+ [(X- V)po + Po(V-X)]VVb 

-povy''>-v[po(A:- v)T/o]} 

(| Ashbv fc Dreitlein|[l975l ). where uniform elastic constants 
were assumed and all quantities Q such as pressure, density, 
gravitational potential V, etc. have been expanded to first 
order in hij, Q = Qo + Q^^\ where Qo denotes the zeroth- 
order quantity and Q'^' the first-order correction. Due to 
spherical symmetry, however, the eigenfunctions of £ are 
of the same type as in equation (|17|) and (|19p . but having 
different coefficient functions a„;(r), 6„;(r) and Cni{r). In 
case of uniform density, po = const., the spheroidal coef- 
fic ient functions a.„i(r) and bni{r) are computed explicitly 
bv lAshbv fc DreitleinI (|l975l ). Since details of the analytical 
expressions of a„;(r), bni{r) and c„;(r) do not enter and af- 
fect the calculations in Section 5 and 6, we simply assume 
eigenfunctions of the type given by equation (|17p and (|19p . 
regardless of whether effects of self-stress have been included 
or not. 



5 INTERNAL VELOCITY FIELD 

In this section, the (complex) internal velocity field v{x,t) of 
the star due to the external gravitational waves is calculated 
from equation (|ll|l . Since the eigenfunctions of £ form a 
complete orthogonal set, the complex velocity field can be 
expanded as 



{x,t) = ^(-iajiv)^jv(t)*jv(a;)e 

JV 



(21) 



where ^jv denote complex time-dependent expansion coef- 
ficients - the amplitudes of the velocity eigenmodes - and 
N = {nlm; Sot T} summarizes all indices of a particular 
eigenfunction. 

The individual amplitudes An of the velocity eigen- 
modes are determined by solving equation pip for each ve- 
locity eigenmode 



* , (r Q 



: c„i(r)iLYir 



(19) 



(22) 
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(cf. equations (|15p and (121(1 ) separately. In the following, 
weak damping is assumed, ryjv <S i^jv for all A'^, where t^at is 
the damping rate. In this limit one can make the approxima- 
tion \dAN/dt\ <^ ljn\An\- Thus, substituting equation p2p 
into equation (|ll|l and applying this approximation yields: 



Of 



_9, ,2 dAjv i^^t 

df ot 
The equation of motion for a particular eigenmode ampli- 
tude An is now obtained by multiplying this equation by 
^^(cc), integrating over the stellar volume V and using the 
normalization condition, equation (|20|) . to give 

(23) 



AAn . — e 

—J- — h (JmAm = 

dt 

Here we have set 



(jjv = ^ y d^s ■ I?(*jv) = 77JV + i5a;jv, 

where 77JV and 5u>m denote the real and imaginary part, re- 
spectively. Introducing the damped eigenfrequency 

UJN = UJN + 5uJN, 

the solution to equation (|23|l after integration by parts with 
respect to time is given by 



AN{t) 



fl + ^ 



dti d xe 
00 Jv 



()7jV+ii<JN)tl 



where we have assumed that the driving term f{x,t) is 
'turned on' at t = —00. Hence, after inserting equation (|12p . 
one immediately finds 



AN{t) 



fl + ^ 



Ati e 



{77jV+i"iv)tl 



(ti) / d^a;p(r)[*^(a;)]'a 



(24) 



Note that thanks to the long-wavelength approximation 
knowledge of the gravitational wave field is not required for 
evaluating the spatial integral, it can rather be computed 
by using properties of the elastic sphere only. The explicit 
computation of the spatial integral and the contraction with 
the gravitational wave tensor is discussed in Appendix A. It 
is found that gravitational waves do solely excite spheroidal 
modes with I = 2, that is ^J[;,„(t) = 0, A^i^{t) = for Z / 2 
and 



4ljJM-l 



dti e 



(77jV+i"N)tl 



hm{ti). (25) 



In order to shorten the notation, the abbreviation = (1 + 
Slon/ujn) was introduced and A*' was set equal to {n2m; S}. 
Furthermore, hm denotes the spherical component of hij 
corresponding to the given azimuthal order m as defined in 
equation ((2]), R denotes the stellar radius and Xn is given 
by the following integral 



dz p{z)z^ [a„2{z) + 3b„2{z)] , 



where z = r/R. Clearly, Xn models the sphere's ability 
to resonate in a particular eigenmode due to the external 
gravitational waves. Note that from equation (|25p . a spe- 
cific eigenmode is only sensitive to the component of the 
gravitational waves having the same azimuthal order m. 



After substitution into equation (|21[l. these amplitudes 
of the velocity eigenfunctions given by equation (f25l) en- 
tirely determine the internal complex velocity field of the 
star under the action of an external gravitational wave field. 
However, the internal (complex) velocity field of a star is not 
a direct observable. As shown in Section 6, relevant quan- 
tities with regard to experiments are the square amplitudes 
|j4jv(f)p of the velocity eigenmodes, rather than the am- 
plitudes themselves. Thus, for the remainder of this section 
we proceed by computing the square amplitudes and start 
with multiplying equation (|25p by its complex conjugate to 
obtain 



\An2m{t)\ 



^^2(Ljr-^ / dfidtae^^^'i+^^-^^^e-'^^'^^-'i^ 



X hm{tl)h*^{t2). 



(26) 



By defining new variables, 

_ tl+t2 _ 

- 7i I = I2 — II, 



to 



equation (1261) can be transformed into 



/ dto e 



2(t-to) 



dTK'^{T)e 



2(t-t„) 



where we have introduced the correlation function 
if*»(r) = ft™ (^fo - C (to + . 

It is now assumed that the gravitational wave field is sta- 
tionary, i.e. the correlation function is independent of to, 
= Km{T). Since we are not interested in the transient 
behaviour of the system after the driving force of the gravi- 
tational waves has been 'turned on' at t = —00, but rather 
in the stationary late-time behaviour, we formally consider 
the limit t — > 00. In this limit the integrals over to and r can 
be performed independently of each other to give the final 
result 



\At 



lim \Af,2m{t)\^ 



SujN_y 

UJN J 



1+^-^1 Pm(t^]v), 



(27) 



where Pmiuj) is the Fourier transform of the correlation func- 
tion Kmir), 



dTKmir)e 



Due to the assumption of a stationary gravitational wave 
field, Km{T) is the auto-correlation function of the spherical 
acceleration hm{t) and thus, Pm{ijj) is the power spectral 
density of hm{t). 

Normally, in the case of weak damping the frequency 
shifts Sojn of the eigenfrequencies due to damping are ex- 
pected to be small. Sun <^ ojn, in which case we have 



\A^ 



32rjMUijfI 



-zPm{ljJN), 



(28) 



where, again, A'' = {n2m; S}. Note the very sensitive R^- 
dependence on the radius of the star. However, when ap- 
plying the normalization condition (|18p this will be reduced 
to an effective Ti'^-dependence. Note in particular the re- 
markable simplicity of this result, equations (|27|) and (|28p . 
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Apart from the aforementioned effective 7?^-dependence and 
numerical pre-factors, the square amplitudes of the veloc- 
ity eigenfunctions are simply obtained by multiplying the 
incoming power spectral density of the spherical accelera- 
tions corresponding to a given azimuthal order m with the 
quantity x^, which, as stated above, describes the resonance 
properties of the star for a given radial order n in case of 
excitation by gravitational waves. The appearance of Pmiio) 
in equation (I27|l and (|28|l can be interpreted as an explicit 
manifestation of the tidal force nature of the excitations. 



Note again the characteristic i?**-dependence (effective B?- 
dependence) on the radius of the star present in the ex- 
pressions for the power Pjv and the surface velocity v% ^ . 
The important relations derived in this section, equa- 



tions (|32|) . (|33|) . (|35|) and (|36|), are entirely ana,logou s 
to the expressions obtained bv | Samadi fc Goupill ll200lll. 
ISamadi. Goupil fc Lebretonl (|200lh and Belkacem et al.l 
('2008I. I2OO9F in case of stochastic excitation by turbulent 
convection. 



6 SURFACE VELOCITY FIELD 

As already stated in Section 5, the quantity accessible to 
experiments is the surface velocity field of a star, rather 
than its internal (complex) velocity field. 

We define the mean-square surface velocity v% g for each 
mode A'' = {n2m; 5*} that is excited by 



(29) 



where ujv(a;,t) is the internal complex velocity field of the 
given velocity eigenmode. 



VN{x,t) = ~iuJNAN{t)'^N{x)e 



(30) 



and () denotes time average. Not e that this defin i tion i s 
equivalent to the one employed by iBelkacem et~aLl l|2009l ). 
Inserting equation (|30|l into the definition given in equation 
(f29ll yields 



(31) 



As shown in Section 5, for a stationary gravitational wave 
field the square amplitudes jyljv(i)P become time indepen- 
dent in the late-time limit t 00, in which case by setting 

*ir(r) = aliir) + Qbliir) 

and inserting equation (|27|l we obtain 



(33) 



VN,a = ^-j p5 IH *]v(«)- (32) 

Again, for small frequency shifts, Sujn <C ujn, this expression 
reads 

Alternatively, expressing the mean-square surface velocity 
in terms of the power that goes into each mode, 

/ 1 
Pjv = 2»7iv( / dm-VN{x,t) ■ VN(x,t) 

= VNO^%l{\AM{t)f), (34) 
one has from equation (I3ip the relation 

^ks = Tph'fliR), (35) 

where in the late-time limit, for small frequency shifts and 
hence with the help of equations (|28|) and (|34p the power is 
explicitly given by 

R'xl 



Pn = 



32/ 



■Pm(t»J]v). 



(36) 



7 CONCLUSIONS 

In the present paper, an excitation mechanism has been es- 
tablished that yields the oscillation amplitudes of stellar os- 
cillation modes when they are excited by arbitrary external 
gravitational waves. This mechanism is expected to be rel- 
evant for solar and solar-like g modes, especially in case of 
stars larger than the Sun that are nearby strong gravita- 
tional wave sources. 

Starting from the general relativistic field equations, the 
inhomogeneous equation of motion for the internal (com- 
plex) velocity field of a non-relativistic star under the ac- 
tion of arbitrary external gravitational waves was derived. 
In this derivation, several assumptions and approximations 
have been made. It was assumed that the linearized theory of 
relativity is applicable and that the internal motions of the 
star are non-relativistic. Furthermore, the long-wavelength 
approximation was employed, which in practice has to be 
verified case by case. Finally, we adopted the isotropic elastic 
sphere model of a star according to the general idea of com- 
bining the physics of spherical resonant mass gravitational 
wave detectors with asteroseismology. This approximation 
was made for calculative simplicity and to demonstrate the 
formalism of the proposed excitation mechanism. 

The equation of motion for the internal velocity field 
was solved analytically and the calculations somewhat 
parallel the methods used in gravitational wave detector 
physics and those in earlier works on excitation mechanisms 
(Balmfort^ 1 19921 : IColdreich eral.i 1994; Sam adi fc Goupill 
200 1|). It is found that external gravitational waves solely 
excite quadrupolar spheroidal normal modes and that eigen- 
modes are sensitive only to the spherical component hm of 
the gravitational waves having the same azimuthal order 
m. The main result of this paper is the expression for the 
stationary square oscillation amplitudes given by equations 
H27|) and (|28|) and was derived under the assumption of a sta- 
tionary external gravitational wave field. The square ampli- 
tudes were shown to have a characteristic _R*-dependence on 
the radius of the star (effective _R^-dependence), which en- 
dows this excitation mechanism with increasing importance 
for stars larger than the Sun, especially for those that are 
nearby strong gravitational wave sources. Furthermore, it is 
worth noting the remarkable simplicity of the main result: 
stationary square amplitudes are obtained by simply mul- 
tiplying a factor that depends on the resonance properties 
of the star and that has an effective i?^-dependene with the 
power spectral density of the spherical accelerations of the 
gravitational waves incident on the star. The latter factor 
can be seen as an explicit manifestation of the excitations' 
tidal force nature. 

Using the result for the square amplitudes, ana- 
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lytical expressions for the mean-square surface veloci- 
ties were derived, which ar e anal og ous to the ones ob - 
tained by Samadi & Goupil ('200l|), ISamadi et aD (j200ll ) 
and ,Belkacem et al., (|2008, 2009) in case of stochastic exci- 
tation by turbulent convection. These surface velocities are 
crucial for assessing the astrophysical and asteroseismic rele- 
vance of the proposed excitation mechanism. It is important 
to note, however, that the mean-square surface amplitudes 
derived in Section 6 are still theoretical quantities. In order 
to predict surface amplitudes due to the present excitation 
mechanism by gravitational waves that can be compared to 
observational data, one has to consider disc-integrated ap- 
parent surface velocities, which take both geometrical and 
limb-darkening effects into account. Since our calculations 
in the present paper are based on the elastic sphere model 
of a star, we do not proceed as far as to give detailed analyt- 
ical expressions for the apparent surface velocities, although 
all important results derived in this paper such as equations 
((27|) . ((28)) . ((32)) and ((33)) also hold in case of a hydrodynam- 
ical computation. 

In future studies, we will concentrate on a hydrodynam- 
ical framework of this excitation mechanism; preliminary re- 
sults show that the aforementioned relations derived in this 
paper remain valid. This would then allow us to explicitly 
calculate apparent surface velocities and thereby to assess 
the asteroseismic relevance of this excitation mechanism. 
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APPENDIX A: DETAILED DISCUSSION OF 
THE SPATIAL GRAVITATIONAL WAVE 
INTEGRAL 

This appendix is concerned with evaluating the spatial in- 
tegral in equation (|24p . which is of the type 



(Al) 



where hij is in TT gauge. Clearly, there are two cases to be 
considered. 



Al Toroidal modes 

Inserting equation p9|) into equation (|A1() yields the expres- 
sion 

hit) I d'..,.,c„,,„.,.,..AK,;,.. 

Jv 

where e,;pq is the Levi-Civita symbol. Integrating by parts 
and using the symmetry properties of the metric tensor hij 
and of the Levi-Civita symbol when contracting one obtains 
the result, that the expression give n in equ ation (IA1() van- 
ishes in case of toroidal modes (see iMaggiore 2008 ) . There- 
fore, toroidal modes do not couple to external gravitational 
waves. 



A2 Spheroidal modes 

For spheroidal modes one has to consider the expression 

hijit) I ArAQ. p{r)r[ani{T)xiXjYi*^ - ib„i(r)eij,qXj,x-jLqY{'^] , 
Jv 
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which is obtained by inserting (|17p into equation (|A1|I . The 
solid angle is denoted by f2. In the following, we consider 
the contributions from the two terms in brackets separately; 
the expressions arising from the first and second term are 
denoted by I and II, respectively. 

By making use of equation ([6]), the first contribution (I) 
can be rewritten as 



In this expression the contribution arising from the first term 
in brackets vanishes, since with the help of equation ([1} one 
has hijSij = h^i = 0. The contribution arising from the sec- 
ond term is easily seen to be non-vanishing only for 1 = 2 
due to the orthogonality property of the spherical harmon- 
ics. Consequently, in order to obtain a simple analytical ex- 
pression for I, we again start with the original expression 
but now setting 1 = 2, 



ArAQ p{r)ra„2{r)hijXtXjY2rr^ 



and substitute equation ^ therein. Then, again using the 
orthogonality property of the spherical harmonics, one is left 
with 

l = hm{t) / Ar p{r)r a„2{r). 
Jo 

Finally we define a new variable, z = r/R, such that the 
integral is independent of the stellar radius R. This yields 

I = h,^{t)R* [ Azp{z)z^a„2iz) = hm{t)R\^^\ 
Jo 

Proceeding as in the former case, we employ equation 
((6]) to gain the following expression for the second contribu- 
tion (II): 



hij / ArAQp{r)r b„i{r)i 
Jv 



borious algebra and again using equations ((3)1 and ([5}: 

2 

II = -2h^R\l^'> + 2j2-h,„,yfp' {yf^y 

m' = -2 

X / drp(r)r^6„2(r) / AQujUp, 



(2) 

where Xn is defined by 



(2) 



Azp{z)z b„2{z). 



With the help of the general identity 
4 



AQ, nifij 



and the orthogonality relation Q the second contribution 
II can finally be written as 

II = 3/i,ni?*Xn'- 

Consequently, it was shown that the expression HA1|I 
is non- vanishing only for spheroidal modes with / = 2, in 
which case 



A'xp{r)[^*j^{x)Yx' 



— hmR Xm 



where we defined 



Xn = Xn +3xi - / Azp{z)z [an2{z) + 'ihn2{z)] . 

Jo 

Note that given a particular density profile p{z) and coeffi- 
cient functions a„2{z) and bn2{z), this quantity is indepen- 
dent of the stellar radius R. 



Since LqYim gives a linear combination of spherical harmon- 
ics with the same value of I, the contribution from the first 
term in brackets vanishes unless I = due to the orthogonal- 
ity property of the spherical harmonics. However, for I = 
the coefficient function 6„; (r) is zero, bno (r) = (see Section 
4). Likewise, the contribution arising from the second term 
in brackets is seen to vanish unless I = 2. 

Knowing that II is non-vanishing only for I = 2, we 
start again with the original expression. 



H] I ArAil p(r)rb„2{r)ei dbY2 

n 

'v 



where Lq = —leqabXadb has already been inserted. Making 
use of the cyclic property of the Levi-Civita symbol and 
contracting the epsilon tensors yields 

hij j ArAQ, p{r)rbn2{r) [r^XjdiY2 dpY2m], (A2) 

Jv 

where the second term in brackets vanishes, since XpdpY2m = 
r-I^Y^^ = 0. Substituting (y^™)*^, (^) for diY^^ in 
equation (|A2|) (cf. equation (O), one obtains after some la- 



